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We present a determination of the CKM matrix element \ V u b \ based on the analysis of semileptonic 
B decays from a sample of 88 million T(4S) decays collected with the BABAR, detector at the PEP- 
II e + e~ storage ring. Charmless semileptonic B decays are selected using measurements of the 
electron energy and the invariant mass squared of the electron-neutrino pair. We obtain |V u i,| = 
(4.41 ± 0.30 to^7 ± 0.28) x 10~ 3 , where the errors represent experimental uncertainties, heavy quark 
parameter uncertainties, and theoretical uncertainties, respectively. 

Note: this document includes corrections to the original publication of this work Q]. These 
corrections, which affect the calculated efficiencies and quantities derived from them, will appear in 
an erratum. 



PACS numbers: 13.20.He, 12.15.Hh, 14.40.Nd 

The study of the weak interactions of quarks has played 
a crucial role in the development of the Standard Model 
(SM), which embodies our understanding of the funda- 
mental interactions. The increasingly precise measure- 
ments of CP asymmetries in B decays allow stringent 
experimental tests of the SM mechanism for CP viola- 
tion via the complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix 0. Improved determinations 
of \V u b\, the coupling strength of the b quark to the u 
quark, will improve the sensitivity of these tests. 

Two observables have been used to determine \V u b\ 
from inclusive semileptonic B decays: the endpoint of the 
lcpton momentum spectrum |3| and the mass of the ac- 
companying hadronic system |4j. In this paper, semilep- 
tonic B — > X u ev decays are selected using a novel ap- 
proach based on simultaneous requirements for the elec- 
tron energy, E e , and the invariant mass squared of the eU 
pair, q 2 [5j. The neutrino 4- momentum is reconstructed 
from the visible 4-momentum and knowledge of the e + e~ 
initial state. The dominant charm background is sup- 
pressed by selecting a region of the q 2 -E e phase space 
where correctly reconstructed B — > X c ev events are kine- 
matically excluded. Background contamination in the 
signal region is due to resolution effects and is evaluated 
in Monte Carlo (MC) simulations. Theoretical calcula- 
tions are applied to the measured B — ► X u eU partial 
rate to determine \V u b\, the precision of which is limited 
mostly by our current knowledge of the &-quark mass, 
m b . 

The data used in this analysis were collected with 
the BABAR detector at the PEP-II asymmetric-energy 
e + e _ storage ring. The data set consists of 88.4 mil- 
lion BB pairs collected at the T(4S) resonance, cor- 
responding to an integrated luminosity of 81.4 fb _1 at 
Js = 10.58 GeV. An additional 9.6 fb _1 of data were 
collected at center-of-mass energies 20 MeV below the 
BB threshold. Off-resonance data are used to subtract 
the non-BB contributions from the data collected at 
the T(4S") resonance. To do so, the off-resonance data 
are scaled according to the integrated luminosity and 
the energy dependence of the QED cross-section, and 



the particles are boosted to the T(45') resonance energy. 
Throughout this paper, all kinematic variables are given 
in the T(4S) rest frame unless stated otherwise. 

The simulation of charmless semileptonic B decays 
used in optimizing the analysis and determining recon- 
struction efficiencies is based on the Heavy Quark Ex- 
pansion (HQE) including O(as) corrections pfl. This 
calculation produces a continuous spectrum of hadronic 
masses, mx- Subsequent hadronization is simulated us- 
ing JETSET down to Im^ |8j. Decays to low-mass 
hadrons (n, r\, p, uj, rj') are simulated separately using 
the ISGW2 model |9j, and mixed with the non-resonant 
states so that the mx, q 2 and E e spectral distributions 
correspond as closely as possible to the HQE calculation. 

Hadronic events containing an identified electron with 
energy 2.1 GeV < E e < 2.8 GeV are selected. Radia- 
tive Bhabha events rejected using the criteria given in 
Ref. and electrons from J/?/> — > e + e~ decays are ve- 
toed. The total visible 4-momentum, p v - ls , is determined 
using charged tracks emanating from the collision point, 
identified pairs of charged tracks from K® — > 7r + 7r~, 
A — > pir~ and 7 — > e + e _ , and energy deposits in the 
electromagnetic calorimeter. Each charged particle is as- 
signed a mass hypothesis based on particle identification 
information. Calorimeter clusters unassociated with a 
charged track and with a lateral energy spread consistent 
with electromagnetic showers are treated as photons. 

Additional requirements are made to improve the qual- 
ity of the neutrino reconstruction and suppress contri- 
butions from e + e~ — > qq continuum events. We form 
the missing 4-momentum, p m - lss — p e + e - — p V is, where 
Pe+e- is the 4-momentum of the initial state. For each 
event we require (1) no additional identified e or /i; (2) 
—0.95 < cos6* m i ss < 0.8, where flmiss is the polar angle of 
the missing 3-momentum; (3) 0.0 GeV < B m i ss — |p miss | < 
0.8 GeV, where B m ; ss is the missing energy in the event; 
(4) |p miss | < 2.5 GeV and (5) | cos6» T | < 0.75, where 9 T is 
the angle between the electron momentum and the thrust 
vector of the remaining particles in the event. 

The measured |p m i ss | differs from the true neutrino 
momentum due to additional particles that escape de- 
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tection. Therefore, a bias correction, p^ = p m i ss (0.804 — 
0.078/|p m i SS |), is derived from the simulation. Since 
the resolution on |p m i ss | is superior to that of i? m i SS , 



we set p v = (p„,lp„|) and q 2 = 
r?± = \/(l ± 0)/ (1 T 0), where /? 



(Pe + Vvf- Defining 
~ 0.06 is the velocity 
of the B meson in the T(4S) frame, the maximum kinc- 
matically allowed hadronic mass squared for a given E e 
and q' 
±2E f 



q 2 - 2m B [E e j]- + q 2 r}+/AE e ) for 



1 r}± , and 



,+q — 2m i 



other- 



no ; no B — > X c ev 



wise. We require s™ ax < 3.5 GeV 
decays can have values of s™ ax below this limit before 
accounting for resolution. The requirements on E e and 
s™ ax and criteria (l)-(5) were chosen to minimize the to- 
tal (experimental and theoretical) expected uncertainty 
<r{\V ub \)/\V ub \. 

The quality of the neutrino reconstruction is eval- 
uated using a control sample (DeV) consisting of the 
decays B — > D°eV(X), where kinematic criteria re- 
sult in the X system typically being no more than a 
7r or 7 from a D* — > D°X transition. The D° is re- 
constructed in the K~tt + decay mode and we require 
\p D o\ > 0.5 GeV and E e > 1.4 GeV. The D° e com- 
bination must satisfy —2.5 < cos#B.£> e < 1.1, where 



cosObdc = (2EbEdc 



*L)/(2|Pfl||Px>e|) is the 



cosine of the angle between the vector momenta of the B 
and the D° e system assuming the only missing particle 
in the B decay was a single neutrino. After the combi- 
natorial background is subtracted using D° mass side- 
bands, the selected sample consists primarily (~ 95%) 
of B — > D°eV and B — > D*eV decays. The control sam- 
ple selection makes no requirements on the other B in 
the event, and can therefore be used to study the impact 
of the modeling of the other B on the neutrino recon- 
struction. Since the unreconstructed X system in the 
B — ► D°eV(X) decays carries away little energy, a good 
estimate (r.m.s.^ 0.2 GeV) of the neutrino energy can be 
obtained from the known B energy and the measured 
D° and e energies, E®°. A second estimate of the neu- 
trino energy is constructed from the visible momentum 
as described previously. Subtracting the first estimate 
from the second gives the distribution shown in Fig. Q 
where the criteria (l)-(5) described above have been im- 
posed. We find good agreement between data and MC; 
the average (r.m.s.) is 0.066 GeV (0.366 GeV) for data 
and 0.072 GeV (0.365 GeV) for simulated events. 

The DeV control sample is also used to improve the 
modeling of the B — > X c eV decays. After relaxing the 
cos 9 B De requirements and subtracting continuum and 
combinatorial backgrounds, we perform a binned x 2 fit to 
the DeV sample in the variables |pd|, E e and cos^ De- 
The fit determines scale factors for the MC components 
B — ► DeV, B — > D*eV and other contributions (85% of 
which are decays to D** states), while keeping the to- 
tal B — > X c eV branching fraction fixed to the measured 
value 11]. The fit increases the B — > DeV and B — > 
D*eV branching fractions to 2.29% and 6.02% (2. 
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FIG. 1: The difference between the two neutrino energy es- 
timates described in the text for continuum-subtracted data 
and simulated BB events for the DeV control sample. 



and 6.52%) for neutral (charged) B mesons, respectively, 
while decreasing the remaining contributions. By design, 
these revised branching fractions respect isospin symme- 
try and are used in the determination of the background. 

Two control samples are used to reduce the sensitivity 
of the efficiency and background estimates to details of 
the simulation: the DeV control sample described above, 
but with E e > 2.0 GeV; and events satisfying the normal 
selection criteria but having s™ ax > 4.25 GeV 2 , a sam- 
ple with < 5% signal decays. Efficiencies e^ at ^ and 
are calculated separately in data and MC as the ratio 
of DeV candidates satisfying criteria (l)-(5) to the total 
DeV sample. The B — > X u eV signal efficiency is multi- 
plied by the ratio of these efficiencies to reduce sensitivity 
to details of the simulation. The s™ ax > 4.25 GeV 2 side- 
band region is used to normalize the simulated s™ ax dis- 
tribution to the data, reducing sensitivity to background 
normalization uncertainties. 

We determine a partial branching fraction 
AB(E,§^) = B(B -> X u eV)f u , unfolded for de- 
tector effects. The acceptance, f u , is the fraction of 
B — > X u eV decays in the region of interest, E e > 2.0 GeV 



and s™ ax < 3.5 GeV , where E e and s™ ax are the true 
(generated) values in the B meson rest frame. Slightly 
lower values are accepted for E e than for E e to account 
for the boost of the B meson and to increase /„. The 
efficiency times acceptance for B — > X u eV decays can be 
written as e u = e S i g /« + e^(l - /„), where e s i g (e^g) 
is the efficiency for an event inside (outside) the region 
of interest to be reconstructed and pass our selection 
criteria. We calculate the partial branching fraction as 
follows: 



2 N 



BB 



1 - fu ' 



fit 



(1) 



where V can d and iV s ide refer to the number of candi- 
dates in the signal and s™ ax sideband regions of the data, 



Mi 



bkg 



and M s i(je refer to background in the signal re- 
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gion and the yield in the sideband region in simulated 
events and 2N BB is the number of B mesons produced 
from T(4S) — > BB decays. Since the resulting ratio of 
e il^/ e sig is small, AB depends only weakly on the model 
used to determine f u . 

Fig.|21shows the electron energy and s™ ax distributions 
after cuts have been applied to all variables except the 
one being displayed. The discrepancy observed between 
data and MC for E e < 1.95 GeV is covered by the sys- 
tematic error on the B — > X c ev modeling. The yields 
and efficiencies are given in Table H] We find 



AS(2.0, 3.5) = (4.41 ± 0.42 ± 0.42) x 10" 



(2) 



where the uncertainties are statistical and systematic, re- 
spectively. Alternative values of AB are obtained using 
different electron energy requirements: A£>(1.9,3.5) = 
(5.29 ± 0.44 ± 0.72) x 10~ 4 and AS(2.1,3.5) = (3.68 ± 
0.43 ±0.36) x 10- 4 . 
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FIG. 2: The electron energy, E e , and s™ ax spectra in the 
T(45) frame for continuum-subtracted data and simulated 
BB events satisfying all selection requirements except for the 
variable shown. The arrows denote the signal (and sideband) 
region in E e and s™ ax . Note that E e / E e (see text). 



Systematic uncertainties are assigned for the modeling 
of the signal B — ► X u ev decays, background and detector 
response. The leading sources of uncertainty are listed in 
Table |H] Uncertainties from the simulation of charged 
particle tracking, neutral reconstruction, charged par- 
ticle identification, and the energy deposition by K\ 
were evaluated from studies comparing data and simula- 
tion. Radiation in the decay process was simulated using 



TABLE I: Yields and efficiencies from data and simulation. 
All uncertainties are statistical except for N B g where system- 
atics are included. Efficiencies are quoted in units of 10 -4 . 



Ac an d 

5130 ± 150 


A sidc 
6152 ± 130 


data 
e Del7 

902 ± 39 


N Bn (I0 b ) 
88.35 ±0.97 


M bkg 
3176 ± 35 


M sidc 
6423 ± 49 


t DeT7 

906 ± 19 


256 ±3 


e sig 

2.7 ±0.2 


fu 

0.174 



PHOTOS [I3 ; comparisons with the analytical result of 
Ref. 0| were used to assess the systematic uncertainty. 
The uncertainty due to bremsstrahlung in the detector 
was evaluated using the method of Ref. ^l| . The uncer- 
tainty in modeling the background was first evaluated by 
varying the total B — > X c ev, B —* Dev and B — > D*ev 
rates within their measured range. Furthermore, the 
form factors for B — > Dev Q and B — > D*ev were 
varied within their uncertainties, and the composition of 
the D** states was modified to include only narrow reso- 
nances, broad resonances, or Goity-Roberts decays (lif : 
the effect of these variations is reduced by the fit to the 
Dev control sample. The modeling of D decays was var- 
ied based on the measurements reported in Ref. ^3] ; the 
variation in the D — > K° L X branching fractions domi- 
nates the uncertainty. The modeling of B — > X u ev de- 
cays is sensitive to the resonance structure at low mass. 
The branching fractions of B — > (-7T, p, uj, 77, rj')ev were 
varied as follows: tt: ±30%; p: ±30%; w: ±40%; simul- 
taneously r] and 77': ±100%. 

TABLE II: Uncertainties on \V ub \ and A/3. 



Source 


<7(|Ki,|)/|Kb| (%) 


a(AB)/AB{%) 


Tracking 


±0.8 


±1.5 


Neutrals 


±1.7 


±3.4 


Electron ID 


±0.5 


±1.0 


Hadron ID 


±1.0 


±2.0 


Bremsstrahlung 


±1.0 


±2.0 


Kl 


±1.3 


±2.6 




±0.6 


±1.1 


Radiation 


±1.9 


±3.8 


B — > A c e77 modeling 


±2.5 


±5.0 


B — > X u eV resonances 


±2.2 


±4.4 


Statistical 


±4.7 


±9.3 


Total experimental 


±6.7 


±13.3 


Heavy quark parameters 
Theoretical 


+14.6 
-10.6 
±6.3 


±1.5 



TB = 



We extract \V wb \ = [AB/{AC, t b )] 1/2 using 
1.604 ± 0.023 ps O. The normalized partial rate, AC 
computed in units of Ar/|U„b| 2 , is taken from Ref. |18j . 
in which the leading terms in the HQE of the B — > 
X u ev spectra are computed at next-to-leading order, and 
power corrections are included at 0(as) for the lead- 
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ing shape function (SF) and at tree level for sublead- 
ing SFs. The values used for the heavy quark parame- 
ters, m b = 4.61 ± 0.08 GeV and /j* = .15 ± 0.07 GeV 2 , 
with a correlation coefficient of —0.4, are based on fits 
to B — > X c iVi moments |l9[, translated to the shape- 
function scheme of Ref . [2(| • 

We find \V ub \ = (4.41 ± 0.30 ±° f 7 ± 0.28) x 10" 3 for 
E e > 2.0 GeV, where the errors represent experimental, 
heavy quark parameters, and theoretical uncertainties, 
respectively. The latter include estimates of the effects 
of subleading SFs [2l|, variations in the matching scales 
used in the calculation, and weak annihilation [22j. No 
uncertainty is assigned for possible quark-hadron duality 
violation. The determination of \V u b\ is limited primarily 
by our knowledge of m b . An approximate dependence 
is \V ub (m b )\ = \V ub (m )\ (1 + 7 (m b ~ m a ) / m a ) , where 
mo = 4.61 GeV. The sensitivity to higher moments of 
the SF is weak; the change in \ V ub \ when varying /i 2 from 
0.03 to 0.35 GeV 2 with m b fixed is 2%, and the impact 
of using alternative SF parameterizations jl^ is < 2%. 
The overall precision on the above result surpasses that 
of Refs. Q and £|, but is comparable to determinations 
of | V ub | which have become available while this paper was 
nearing completion |24|. 
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